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LAMINAR  BOUNDARY  LAYER  ON  A  CONE  IN 
SUPERSONIC  FLOW  WITH  UNIFORM  MASS  TRANSFER 


Paul  A.  Lfibby 


SUMMARY 

A  solution  for  the  laminar  boundary  layer  on  a  cone  with  uni¬ 
form  mass  transfer  is  obtained.  The  velocity  field  is  found  for  either  suction 
or  injection  but  the  related  solution  for  the  energy  field  is  subject  to  an  energy 
balance  at  the  exposed  surface  and  is  therefore  valid  only  for  injection.  This 
latter  solution  is  equally  applicable  to  certain  species  fields  as  well.  The 
present  results  along  with  those  presented  previously  for  the  two-dimensional 
case  permit  a  comparison  of  the  effect  of  injection  on  boundary  layers  over 
two-dimensional  and  conical  surfaces, 

INTRODUCTION 

The  characteristics  of  laminar  boundary  layers  involving  mass 
transfer  are  of  interest  in  connection  with  the  thermal  protection  of  surfaces 
subjected  to  high-energy  gas  flows  and  with  boundary  layer  control.  The  es¬ 
sential  simplification  resulting  from  the  assumption  of  flow  similarity  has 
led  to  its  widespread  application  to  the  theoretical  treatment  of  such  boundary 
layers.  Indeed,  in  certain  applied  problems,  e.  g.  ,  at  stagnation  points  and 
on  surfaces  involving  steady-state  sublimation,  the  similarity  assumption  is 
exactly  applicable.  Much  less  attention  has  been  devoted  to  nonsimilar  bound¬ 
ary  layers  with  mass  transfer  although  they  arise  in  problems  of  practical 
interest;  e.  g.  ,  experimental  research  of  a  fundamental  nature  is  often  carried 
out  on  porous  surfaces  which  lead  to  uniform  rates  of  suction  or  injection  and 
thus  to  nonsimilar  flows.  The  present  report  deals  with  the  boundary  layer 
on  a  cone  in  supersonic  flow  with  uniform  mass  transfer. 

"^This  report  is  an  extended  version  of  a  paper  ac¬ 
cepted  for  publication  in  the  Physics  of  Fluids.  In  the  abbreviated  version 
the  previous  analysis  of  reference  1  is  relied  on  heavily  for  notation  and  de¬ 
tails;  the  present  extended  version  is  more  complete  and  self-contained. 


The  analysis  develops  first  the  velocity  field  in  the  boundary 
layer  on  the  cone;  this  will  be  found  to  be  applicable  to  suction  as  well  as  to 
injection.  However,  the  next  stage  of  the  analysis  pertains  to  the  energy 
and  species  conservation  fields  associated  with  this  velocity  distribution 
under  steady  flow  conditions;  thus  energy  and  species  balance  conditions  are 
applied  at  the  porous  surface  and  the  resultant  solutions  are  applicable  only 
to  the  case  of  injection.  This  is  generally  the  case  of  greater  applied  interest. 

The  present  analysis  should  be  considered  a  generalization  to 
the  case  of  a  cone  of  previous  two-dimensional  solutions  with  uniform  mass 
transfer.  In  the  context  of  modern  transformations  of  the  equations  describ¬ 
ing  laminar  boundary  layers  for  both  two-dimensional  and  axisymmetric  flows, 
it  may  be  a  priori  surprising  that  these  two  cases  must  be  treated  separately. 
However,  the  requirement  of  uniform  mass  transfer  leads  mathematically  to 
two  distinct  cases  and  the  analysis  for  one  cannot  be  interpreted  so  as  to  apply 
to  the  other.  Physically,  the  distinct  behavior  in  the  two  cases  appears  to  be 
related  to  the  different  effect  of  the  uniform  mass  transfer  on  the  boundary 
layer  growth  which  in  both  cases  is  nonparabolic,  a  symptom  of  nonsimilarity. 

The  problem  of  the  two-dimensional  boundary  layer  with  uniform 
mass  transfer  has  been  considered  by  several  authors;  Iglish  (reference  Z) 
presents  a  solution  for  incompressible  flow  with  uniform  suction.  Hew  and 
Fannucci  (reference  3)  provide  the  extension  to  compressible  flows  with  either 
uniform  suction  or  injection.  The  analyses  of  both  references  Z  and  3  deal 
primarily  with  the  velocity  field  and  involve  the  numerical  solution  of  a  par¬ 
tial  differential  equation  which  is  free  from  parameters  and  which  thus  has 
been  integrated  once-for-all. 

Libby  and  Chen  (reference  4)  recently  treated  the  two-dimensional 
case  by  carrying  out  an  expansion  in  terms  of  a  mass  transfer  parameter. 

This  permits  the  velocity,  energy  and  species  fields,  the  latter  two  involving 
balance  conditions  at  the  porous  surface,  to  be  determined  from  the  successive 
solution  of  ordinary  differential  equations.  The  velocity  field  provided  by  ref¬ 
erence  4  duplicates  that  obtained  previously  in  references  Z  and  3  but  its  avail¬ 
ability  and  describing  form  permit  the  energy  and  species  fields,  which  had 
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previously  been  l«±rgely  ignored,  to  be  re3.dily  obtained.  Of  interest  in  these 
solutions  i$  the  distribution  of  wall  enthalpy  and  wall  concentration  of  the 
various  species  present.  For  example,  at  the  leading  edge  of  the  plate  the 
enthalpy  at  the  porous  surface  is  equal  to  the  stagnation  enthalpy  of  the  ex¬ 
ternal  flow;  downstream  from  the  leading  edge  the  wall  enthalpy  approaches 
that  in  the  coolant  chamber.  Similar  considerations  apply  to  the  composition 
at  the  porous  surface.  Such  streamwise  behavior  is,  of  course,  symptomatic 
of  nonsimilar  flows.  No  generalization  of  the  aforementioned  two-dimensional 
analyses  to  the  axisymmetric  case  appears  to  have  been  performed.  However, 
Lribby  (reference  4)  in  connection  with  an  analysis  of  mass  transfer  effects 
treated  as  perturbations  recently  performed  a  calculation  related  to  that  pre¬ 
sented  here  but  involving  only  the  first-order  effect  of  injection.  The  prin¬ 
cipal  motivation  and  utility  of  the  analysis  of  reference  4  are  the  possibility 
of  computing  the  effects  of  arbitrarily  distributed  mass  transfer. 

It  is  perhaps  of  interest  to  note  that  none  of  the  aforementioned 
analyses  are  able  to  answer  the  question  of  whether  a  constant,  finite  injection 
rate,  (pv)  ,  applied  over  a  finite  streamwise  length,  x,  leads  to  ’’blow-off,  ’’ 
i.  e.  ,  to  a  point  of  zero  skin-friction.  The  earlier  studies  of  boundary  layers 
with  mass  transfer  based  on  momentum  integral  methods  indicated  a  negative 
answer  to  this  question  but  Stewartson  (reference  5)  has  recently  called  at- 
tention  to  some  unpublished  results  which  imply  an  affirmative  answer. 

Studies  of  the  behavior  of  boundary  layers  in  the  neighborhood  and  downstream 
of  the  "blow-off”  point  are  of  considerable  interest  but  are  outside  the  scope 
of  the  present  paper. 

The  present  analysis  follows  closely  that  of  reference  1  and 
provides  its  extension  to  the  axisymmetric  case;  the  solution  for  the  velocity 
field  is  presented  first  and  then  applied  to  the  determination  of  the  energy 
and  species  fields. 


‘in  terms  of  the  mass  transfer  parameter  of  reference  1,  namely 

€  =  [“(pv)  /p  p  u  ](s/2)^^^,  "blow-off"  occurs  at  e  =  -0.606. 

,  w  e  e  e 
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SYMBOLS 


r 

s 

u 

V 
X 

Y. 

1 

y 

Of 

c 

V 
Qc 

p 

X 

Subscripts 


local  skin-friction  coefficient 

transfornaed  stream  function,  cf.  equation  (Z) 

unit  solutions  for  the  energy  and  species  fields,  cf.  equation  (ZO) 

ratio  of  stagnation  enthalpies,  h  /h 

s  s ,  e 

static  enthalpy 

stagnation  enthalpy,  (u^/Z)  +  h 

unit  solutions  for  the  velocity  field,  cf.  equation  (11) 
inhomogeneous  terms,  cf.  equations  (IZ)  and  (Zl) 
cylindrical  coordinate 

transformed  streamwise  coordinate,  cf.  equation  (1) 

streamwise  velocity  component 

normal  velocity  component 

streamwise  coordinate  along  the  surface 

mass  fraction  of  species  i 

coordinate  normal  to  the  surface 

cone  parameter,  Oi  =  sin  0 

c 

mass  transfer  parameter  for  two-dimensional  case 
transformed  normal  coordinate,  cf.  equation  (1) 
cone  half-angle 
viscosity  coefficient 
mass  density 

transformed  streamwise  coordinate,  cf.  equation  (8). 


c  refers  to  conditions  in  the  coolant  chamber 

e  refers  to  conditions  in  the  external  stream 

w  refers  to  conditions  at  the  surface,  y  ^  rf  =  0 

^  refers  to  stagnation  conditions  in  the  external  stream. 


THE  VELOCITY  FIELD 


The  geometrical  and  dynamical  variables  being  considered  are 
shown  schematically  in  figure  1.  The  velocity  distribution  throughout  the 
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boundary  layer  may  be  obtained  in  terms  of  transformed  variables  from  the 
x-wise  momentum  and  the  mass  conservation  equations  alone  provided  there 
is  employed  the  common  assumption  that  the  ratio  1.  The  trans¬ 

formed  variables  are  the  well-known  Levy-Lees  Tj,  s  (cf.,  e.  g.  ,  references 
5  and  6)  defined  for  the  present  problem  according  to 


^  =  Pe  ^  Jq 


2  j  2  ,  3,,, 

s=puu  ;  r  dx=p/iu  a  (x/3) 


e'  e  e 


(1) 


where  r  =  x  sin  6  =  ax.  The  velocity  components  are  found  from  the  trans- 

c 

formed  stream  function  f(s,  tj)  by  means  of  the  equations 


e  TJ 


(^) 


=  -  [(2s)^^^/pr]ip^PgU^r'^  [(f/2s)  +  f^]  +  f^}  . 


A  useful  result  from  the  second  of  equations  (2)  relates  the  distributions  of 

(pv) 


f(s,0)  =  f  (s)  and  of  (pv)  ,  i.  e.  , 
w  w 


d  r,,  .1/2  f  W 

ds  w^  ■  pUur 


(3) 


e'  e  e 


Finally,  the  partial  differential  equation  for  f(s,T))  is  (cf.  ,  e.  g.  ,  reference  6) 


rjTJT? 


+  f  f 


7]T? 


<V^S77'/s  V 


(4) 


Consider  next  the  conditions  to  which  the  solution  of  equation  (4) 

is  to  be  subject;  if  it  is  assumed  that  s  f  0  for  all  TJ  b.s  s  ^  0  and  that 

s  j  y  ^ 

f(0,  0)  =  f^(0)  =  0  as  will  be  shown  to  be  the  case  below,  indeed  f^  s  , 
then 


yo.rj)  =  f^(v) 


(5) 
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where  is  the  Blasius  function  and  where  the  prime  denotes  differentiation 
with  respect  to  77.  This  behavior  near  the  origin  of  the  boundary  layer  should 
be  contrasted  with  that  which  prevails  for  flows  with  f^{0)  ^  0;  this  is  the 

case  for  the  boundary  layer  on  a  cone  in  supersonic  flow  if  (nv) 

^  w 

at  least  as  s  0.  Then  the  effect  of  mass  transfer  prevails  even  at  s  =  0. 

In  the  present  case  the  uniform  mass  transfer  must  persist  over  some  dis¬ 
tance  from  the  origin  before  it  influences  the  boundary  layer. 

The  conditions  on  f^{s,rj)  far  from  and  at  the  porous  surface 

are  the  usual 


f^{s,cx>)  =  l  ,  f^(s,0)=0  .  (6) 

The  final  condition  introduces  the  effect  of  mass  transfer;  from  equation  (3) 

it  is  found  that  for  uniform  mass  transfer,  i.  e.  ,  (pv)  =  constant. 

'w  ’ 

f(s.O)  =  f^{s)  =  {L-(pv)^/p^u^]  . 

(7) 

®  contrasted  to  the  two-dimensional  case  in  which 

1  /  Z 

®  to  the  case  of  similar  flow  in  which  f  (s)  =  constant.  For 

w 

a  given  uniform,  injection  rate,  given  cone  angle  and  given  properties  of  the 
external  flow  the  factor  within  {  ]  in  equation  (7)  is  known. 


It  is  convenient  to  introduce  a  new,  nondimensional,  independent 
variable,  such  that  the  solution  for  f  is  considered  as  f{x»  V)  where 

X  =  {  }  .  (8) 

Note  that  X  <  0  for  injection  and  X  >  0  for  suction.  Note  further  that  x 
corresponds  to  e(s)  used  in  reference  1  for  the  two-dimensional  case;  this 
correspondence  will  be  of  interest  in  some  considerations  below  of  the  relative 
effect  of  mass  transfer  on  the  boundary  layer  on  wedges  and  cones. 

Equations  (4)  through  (7)  then  yield  the  following  mathematical 

problem: 
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f  +  f  f  -  (v/3)  (f  f  -  f  f  )  =  0 

mv  m  ^  V  nx  x  m 


yO.T?)  =  VqW  :  yx.oo)  =  1  ;  yx.o)  =  0  ;  f(x,  0)  =  x  • 


(^) 


The  quantity  of  most  technological  interest  from  the  solution 
for  the  velocity  field  is  the  skin-friction  coefficient;  in  terms  of  f(x.  Tj)  this 
is 

=  (10) 

SO  that  the  distribution  of  wall  values  of  f _ is  of  interest.  Note,  of  course, 

Tjr? 

that  in  equation  (10)  X  “  X(^)* 

The  solution  of  equations  (9)  can  be  carried  out  by  finite  dif¬ 
ference  procedures  for  both  positive  and  negative  values  of  X  starting  at 
X  =  0.  Such  a  solution  would  correspond  to  the  extension  to  the  cone  case 
of  the  two-dimensional  calculations  of  Iglish  (reference  2)  for  X  ^  ^ 
of  Liew  and  Fannucci  (reference  3)  for  X  $  less  ambitious  computing 

task  is  involved  if  a  series  solution  analogous  to  that  of  reference  1  is  found; 
such  a  solution  appears  to  be  adequate  for  many  applied  problems.  Thus 
assume 


oo 

HX’ v)  =  +  Y. 

n=l 

Substitution  into  equation  (9)  and  collection  of  like -powers  in  x  leads  to  the 
infinite  array  of  equations 


N»'  +  ^0  ^n  ^  ^0  ^  Ll  +  (n/3)]  =  0  ,  n  =  1 


n  >  2 


(1^) 


where  the  are  nonlinear,  inhomogeneous  terms  depending  on  N,  (r)), 

n 

k  =  l,  2,  n-1.  For  example, 
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(13) 


-  (1/3)  t(H' )^  *  Nj  N'-  ] 

-  N'^  N^  -  (1/3)  N^  N^  -  (2/3)  N^  N'^  . 

The  solutions  of  equations  (IZ)  are  subject  to  the  boundary  conditions 

n;^(o)  -  n;^(oo)  =  0 

=  1  ,  n  =  1  (14) 

=  0  ,  n  >  Z 

The  first  five  functions  N  (tf)  have  been  found  numerically. 

The  crucial  values  for  generating  these  functions  in  detail  are  N”{0);  ac¬ 
cordingly,  these  have  been  listed  in  Table  I,  In  addition  the  ’’velocity  pro¬ 
files  are  usually  of  most  graphic  interest  so  these  have  been  given  in  fig- 

Also  of  interest  is  the  variation  of  (f  )  with  y:  this  is  shown  in 

r)T7  w 

figure  3  for  four  and  five  terms  in  the  series.  From  a  practical  viewpoint 
the  availability  of  only  five  terms  would  appear  from  figure  3  to  be  unimpor¬ 
tant  provided  |xl  ^  0.  6. 

The  availability  of  this  solution  for  the  velocity  field  in  the  axi- 
symmetric  case  along  with  the  previously  available  solutions  for  the  two- 
dimensional  case  permits  a  comparison  of  the  relative  effects  of  mass  trans¬ 
fer  on  cones  and  wedges  be  carried  out.  For  this  purpose  it  is  of  in¬ 
terest  to  compare  c(s)  and  X(^)i  thus  consider  a  wedge  and  a  cone 
under  flow  conditions  so  that  u^  are  the  same  in  the  two 

cases  and  let  (pv)  /p  u  be  common.  Then  at  equal  distances  from  the 
leading  edge  and  from  the  apesi,  d.  e.  ,  at  the  same  station  x,  it  can  be  shown 
that 

e  =  (2/3^/^)  X  .  (15) 

Moreover,  the  ratio  of  the  skin-friction  coefficients  at  the  same  x-station 
becomes 
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^*'£^cone 


f  (s.O) 

TJZ 


(16) 


where  it  is  to  be  understood  that  the  s  and  x  the  arguments  of  are 

related  according  to  equation  (15)  in  order  to  assure  comparison  at  the  same 

x-station.  Note  that  if  f  (s,  0)  =  f  ^(X>0)»  then  equation  (16)  is  simply  the 

r]7]VV 

well-known  relation  between  the  skin-friction  on  a  cone  and  wedge  without 
mass  transfer.  If  f  {s,0)  =  f  {Yy  0)  with  mass  transfer,  then  on  the  basis 
of  the  comparison  as  carried  out  here,  mass  transfer  is  considered  equally 
effective  in  both  geometries. 


In  figure  4  the  ratio  of  the  two  skin-friction  coefficients  is  given 
for  a  range  of  X-  wide  excursion  of  this  ratio  from  the  zero  mass  transfer 

value  of  0,  578  will  be  noted.  From  this  figure  it  will  be  seen  that  the  effect 
of  mass  transfer  in  altering  the  skin  friction  is  greater  on  wedges  than  on 
cones  when  the  comparison  is  carried  out  on  the  basis  described  above. 

Note,  of  course,  that  the  total  mass  added  through  the  porous  surface  from 
the  apex  of  the  cone  to  the  generic  station  x  per  unit  length  of  perimeter  at 
X  is  one-half  that  added  from  the  leading  edge  of  the  wedge  over  the  same  x- 
wise  length  per  unit  length  parallel  to  the  leading  edge. 

THE  ENERGY  FIELD 


The  energy  field  will  be  described  in  terms  of  the  ratio  of  stag¬ 
nation  enthalpies,  g  =  h  /h  .  In  addition  to  the  previously  employed  as- 

s  s ,  e 

sumption,  ==*  1 ,  it  will  be  assumed  here,  as  is  frequently  done,  that 

the  Prandtl  number  is  unity,  that  a  single  diffusion  coefficient  exists  and  that 
the  Schmidt  number  based  thereon  is  unity.  Then  the  energy  equation  in  x»  V 
variables  is  (cf.  ,  e.  g.  ,  reference  6) 


g  +  f  g  “  (x/3)(f  g  -  f  g  )  =  0  (17) 

rj  X  X  n 

where  it  is  assumed  that  f(x»T3)  is  known  from  the  above  analysis.  The  solu¬ 
tions  to  equation  (17)  apply  to  a  variety  of  physical  cases  depending  on  the  ini¬ 
tial  and  boundary  conditions  at  the  porous  surface.  Here  there  will  be  con¬ 
sidered  only  the  case  of  injection  (x  <  0);  moreover,  it  will  be  assumed  that 
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the  convective  heat  transfer  from  the  gas  to  the  porous  surface  is  absorbed 
by  the  injected  gas  in  passing  from  an  injection  chamber,  where  its  enthalpy 
ratio  is  constant  and  denoted  by  g^,  to  the  exposed  surface  where  its  en¬ 
thalpy  ratio  is  variable  and  denoted  by  g^  =  to  be  determined.  Thus 

under  steady  flow  conditions  the  convective  heat  transfer  is  found  to  be  re- 
lated  to  the  enthalpy  difference  (g^  ”  according  to"^ 

{g„)^  =  {4/3)(-x)(g^  -  g^)  .  (18) 

Note  that  both  the  case  of  a  coolant  being  injected  to  provide  thermal  insula¬ 
tion,  i.  e.  ,  <  1,  and  the  case  of  a  heated  gas  being  injected  to  heat  the 

surface  and  the  boundary  layer,  i.  e.  ,  >  I,  may  be  treated  by  the  present 

analysis;  perhaps  the  former  case  is  of  greater  current  interest.  It  is  perhaps 
worth  noting  that  if  in  addition  to  the  convective  heat  load  there  exists  a  uni¬ 
form  additional  thermal  load,  e.  g.  ,  due  to  radiation,  then  this  heat  balance 
condition  still  prevails  but  with  a  redefined  g^  which  accounts  for  this  added 
load. 

F rom  equation  (18)  and  quite  separately  on  physical  grounds 
it  is  expected  that  at  the  origin  (x  =  0),  g^(0)  =  1  so  that  the  initial  and  ex¬ 
ternal  stream  conditions  are 

g(X»«>)  =  g(0,T?)  =  1  .  (19) 

The  solution  of  the  problem  posed  by  equations  (17)  through  (19) 
again  may  be  found  by  finite  difference  calculations;  here  a  series  solution 
is  found  in  the  form 

oo 

g(x.»7)  =  1  +  (g^  -  1)  ^  x*"  GJn)  (20) 

n  =  l 

where  the  G^(T?)  functions  are  given  by  an  array  of  ordinary  differential 
equations 


Reference  1  provides  a  detailed  derivation  of  the  energy  and  species 
balance  conditions;  little  change  is  necessary  to  obtain  the  conditions  for  the 
cone  case. 
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and  are  subject  to  the  conditions 


G^(0)  -  -(4/3)  G^_^{0)  , 


-  (4/3)  ,  n  =  1  . 

Again  the  functions  are  known,  nonlinear  functions  of  the  previous 

N  and  G  functions;  for  example 
n  n 


-  -  (4/3)  G'^  -  (1/3)  N'j  Gj 


(22) 


=  .  (4/3)  G^  +  (2/3)  N'j  G^  +  (1/3)  Gj  -  (5/3)  G'^  . 

The  first  five  G  solutions  have  been  obtained  numerically; 
n 

their  reproduction  by  straightforward  numerical  integration  requires  the 

values  of  G  (0)  and  G'  (0)  listed  in  Table  II.  For  graphic  displace  the 
n  n 

distributions  of  ^^(rj)  are  shown  in  figure  5  while  the  distribution  of  wall 
enthalpy  in  the  form  (1  -  g^)(^  ■  g^)  ^our  and  five  terms  in 

the  series  is  shown  in  figure  6. 


The  quantity  which  is  of  greatest  applied  interest  from  the  energy 
solution  obtained  here  is  the  distribution  of  wall  enthalpy  in  terms  of  g^(x) 
obtainable  from  figure  6.  Note  that  increasing  |x|  corresponds  to  increas¬ 
ing  X  and  that  as  x  increases  g  changes  from  unity  to  g^  in  the  manner 
shown  in  figure  6.  Again  such  behavior  is  symptomatic  of  nonsimilar  flows. 

If  there  is  considered  the  comparative  efficacy  of  injection  on 
wedges  and  cones  in  making  g  approach  g  ,  then  on  the  same  basis  as 
used  above  for  the  comparison  of  skin  friction,  it  is  found  that  injection  is 
more  effective  on  wedges  than  on  cones. 

In  closing  this  section  several  remarks  are  perhaps  in  order. 

The  solutions  for  both  the  velocity  and  energy  fields  have  been  obtained  in 


terms  of  the  transformed  variables  \  and  77;  the  relation  between  the  phys¬ 
ical,  streamwise  coordinate  x  and  \  is  simply  given  in  terms  of  the  in¬ 
jection  rate  and  of  the  external  flow  and  cone  characteristics.  The  trans¬ 
formation  inverse  to  that  of  equation  (1)  yields  y  =  y(X»  terms  of  an 

integral  with  respect  to  Tj  of  the  density  ratio  p^/o  and  requires  explicit 
determination  of  that  density  ratio.  This  in  turn  requires  simultaneous  con¬ 
sideration  of  the  velocity,  energy,  and  species  fields  and  involves  an  equation 
of  state,  and  relations  between  static  temperature  and  static  enthalpy  of  each 
species  present.  Calculations  of  this  sort  are  straightforward  and  are  of  in¬ 
terest  in  connection  with  comparisons  between  experiment  and  analysis. 

Note  that  the  applicability  of  the  present  solution  for  the  energy  field  to  a 
species  field  is  discussed  below. 

CONCLUDING  REMARKS 


In  reference  1  it  is  shown  that  the  solutions  for  species  and 
energy  conservation  in  the  two-dimensional  case  are  identical  provided  no 
chemical  reaction  takes  place,  and  provided  a  species  mass  balance  across 
the  porous  surface  is  asserted.  The  same  considerations  apply  to  the  cone 
case  as  well.  The  mathematical  formulation  for  the  species  field  involves 
replacing  g  in  equation  (17)  by  the  species  mass  fraction,  Y.,  replacing 
^w’  ^c  energy  balance  of  equation  (18)  by  Y.  ,  Y.  and 

[(Y.)^]^,  respectively,  and  replacing  equation  (19)  by  Y.(x,  00)  =  Y.(0,  t?)  = 
Then  the  solution  for  Y^  is  sought  in  the  form 


Y.(x,r?) 


Y.  +  (Y. 

1,  e  1,  c 


00 


Y. 

1. 


n  =  l 


X""  G^(rj) 


as  in  equation  (ZO).  Thus  the  solutions  for  G^(rj)  may  be  applied  in  deter¬ 
mining  the  species  fields  as  well. 


It  is  also  noted  in  reference  1  and  repeated  here  for  complete¬ 
ness  that  the  analysis  of  species  conservation  is  applicable  to  boundary  layers 
with  chemical  reaction  provided  the  species  conservation  is  replaced  by  ele¬ 
ment  conservation;  in  this  case  the  species  and  temperature  fields  must  be 
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found  from  the  solutions  for  element  and  energy  conservation  obtained  here 
and  from  additional  assumptions  regarding  the  chemical  behavior  of  the  sys¬ 
tem,  i.  e.  ,  whether  it  corresponds  to  equilibrium  chemistry  or  to  finite  rate 
behavior,  and  regarding  the  relations  between  static  temperature  and  species 
enthalpy . 
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Table  I.  Initial  values  of  the  N^(77)  functions. 


n 

N"(0) 

n 

1 

0. 9039 

Z 

0.2504 

3 

-0. 1068 

4 

0. 04038 

5 

-0. 07963 

Table  II. 

Initial  values  of  the 

0^(77)  functions. 

n 

g;(0) 

G  (0) 

1 

1.  333 

-2. 465 

Z 

3.  287 

-1. 645 

3 

2.  193 

-0. 2046 

4 

0. 2728 

-0. 1205 

5 

0. 1607 

-0. 1214 
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Figure  4.  The  relative  effect  of  mass  transfer  on  skin  friction. 
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